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Abstract. The Rayleigh conjecture on the representation of the scattered 
field in the exterior of an obstacle D is widely used in applications. How- 
ever this conjecture is false for some obstacles. AGR introduced the Modified 
Rayleigh Conjecture (MRC), and in this paper we present successful numerical 
algorithms based on the MRC for various 2D and 3D obstacle scattering prob- 
lems. The 3D obstacles include a cube and an ellipsoid. The MRC method is 
easy to implement for both simple and complex geometries. It is shown to be 
a viable alternative for other obstacle scattering methods. 



Key words and phrases: obstacle scattering, modified Rayleigh conjecture, nu- 
merical solution of obstacle scattering problem. 

1. Introduction. 

In this paper we present a novel numerical method for Direct Obstacle Scattering 
Problems based on the Modified Rayleigh Conjecture (MRC). The basic theoretical 
foundation of the method was developed in [7]. The MRC has the appeal of an 
easy implementation for obstacles of complicated geometry, e.g. having edges and 
corners. In our numerical experiments the method has shown itself to be a compet- 
itive alternative to the BIEM (boundary integral equations method), see [8]. Also, 
unlike the BIEM, one can apply the algorithm to different obstacles with very little 
additional effort. 

However, a straightforward numerical implementation of the MRC, may be in- 
efhcient or fail. Nevertheless, one can make the numerical implementation of the 
method based on MRC work by introducing some new ideas within the same the- 
oretical framework. In this paper we describe such an idea (Random Multi-point 
MRC) that made possible, for the first time, to apply the method successfully to 
3D obstacles including a sphere and a cube. 

In our previous paper [8] we described another implementation of the MRC. 
That method (Multi-point MRC) could be used for 2D obstacles of a relatively 
simple geometry, but it failed for some 2D obstacles, and it was not successful for 
3D problems because of the problem size limitations. In this paper, in addition 
to treating 3D problems, we review our earlier results and show that they can be 
significantly improved by using the new idea for the implementaion of the MRC 
algorithm. 

We formulate the obstacle scattering problem in a 3D setting with the Dirichlet 
boundary condition, but the method discussed can also be used for the Neumann 
and Robin boundary conditions. 
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Consider a bounded domain D (ZM.^ , with a boundary S which is assumed to be 
Lipschitz continuous. Denote the exterior domain by D' = M.^\D. Let a, a' € S"^ 
be unit vectors, where S"^ is the unit sphere in M^. 

The acoustic wave scattering problem by a soft obstacle D consists in finding 
the (unique) solution to the problem (1.1)-(1.2): 

(1.1) (V^ + fc^) M = in £)', u = OonS, 

gikr / \ ^ 

(1.2) u = uo + A(a',a) hoi-), r := |a;| ^ oo, a' := —. 

r \r J r 

Here uq := e'*^"'^ is the incident field, v := u — uq \s the scattered field, A{a' ,a) 
is called the scattering amplitude, its k-dependence is not shown, fc > is the 
wavenumber. Denote 

(1.3) Ae{a):= [ A{a' ,a)Ye{^da' , 

where Yi{a) are the orthonormal spherical harmonics, = Yim, —£< m < £. 

Informally, the Random Multi-point MRC algorithm can be described as follows. 

Fix a J > 0. Let a;^, j = 1, 2, J be a batch of points randomly chosen inside 
the obstacle D. For x € D' , let 

(1.4) a' = Mx,Xj) = Yi{a')ht{k\x - xj\), 

\X Xj\ 

where hi{r) are the spherical Hankel functions, normalized so that hi{r) ~ as 
r +00. 

Let g{s) = uo{s), s S S, and minimize the discrepancy 

J L 

(1-5) ^(c) = \\g{s) + J2J2'^^^M^'^j^\\LHS)^ 

j=l l=Q 

over c G C"^ , where c = {cf j}. That is, the total field u{s) = g{s) + v{s) is 
desired to be as close to zero as possible at the boundary S, to satisfy the required 
condition for soft scattering. If the resulting residual r""" ~ min $ is smaller than 
the prescribed tolerance e, then the procedure is finished, and the sought scattered 
field is 

J L 

M^) = '^'^ce,jtpe{x,Xj), xeD', 

j=i e=a 

(see Lemma 2.2 below, and the Remark following it). 

If, on the other hand, the residual r™'" > e, then we continue by trying to 
improve on the already obtained fit in (1.5). Adjust the field on the boimdary by 
letting g{s) := g{s) + v^{s), s G S. Create another batch of J points randomly 
chosen in the interior of D, and minimize (1.5) with this new g{s). Continue with 
the iterations until the required tolerance e on the boundary S is attained. In 
each iteration accumulate new interior points Xj and the corresponding best fit 
coefficients cej . After the desired tolerance is reached, the sought scattered field 
is computed anywhere in D' . 

Here is the precise description of the algorithm. 

Random Multi-point MRC. 



MRC METHOD FOR OBSTACLE SCATTERING PROBLEMS 



3 



For Xj e D, and £>0 functions tpe{x,Xj) are defined as in (1.4). 

(1) Initialization. Fix e > 0, i > 0, J > 0, N^ax > 0. Let n = 0, and 

gi-s) = uo{s), s G S. 

(2) Iteration. 

(a) Let n := n+1. Randomly choose J points xj"^ € £>, j = 1,2, . . . , J. 

(b) Minimize 

$(c) = \\g{s) + Y,j2''^,jMs,4'^)\\L^iS) 
over c e C^, where c = {c^j}. 

Let the minimum of ^ be attained at c'"^ = {c^"'')}, j = 1,2, . . . , J, 
and the minimal value of $ be r™*". 

(3) Stopping criterion. 

(a) If r™™ < e, then stop. Compute the approximate scattered field 
anywhere in D' by 

(1.6) v,{x) := EEE^S'/^'.^^'^f )' ^ e D'. 

k=i j=i e=o 

(b) If r™'" > e, and n ^ N^a:,, let 

,/ L 

and repeat the iterative step (2). 

(c) If r"*" > e, and n = Nmax, then the procedure failed. 



2. Direct scattering problems and the Rayleigh conjecture. 

Let a ball Br ■= {x : \x\ < R} contain the obstacle D. In the region r > R the 
solution to (1.1)-(1.2) is: 



X 

? 

r 



(2.1) u{x,a)=e''"^-'' + Y,M(^)^t, ■.= Yi{a')hi{kr), r > R, a' 

where the sum includes the summation with respect to w, —£< m < i, and Ai{a) 
are defined in (1.3). 

The Rayleigh conjecf,ure (RC) is: the series (2.1) converges up to the boundary 
S (originally RC dealt with periodic structures, gratings). This conjecture is false 
for many obstacles, but is true for some ([1, 3, 9]). For example, if n = 2 and D is 
an ellipse, then the series analogous to (2.1) converges in the region r > a, where 
2a is the distance between the foci of the ellipse [1]. In the engineering literature 
there are numerical algorithms, based on the Rayleigh conjecture. Our aim is to 
give a formulation of a Modified Rayleigh Conjecture (MRC) which holds for any 
Lipschitz obstacle and can be used in numerical solution of the direct and inverse 
scattering problems. We discuss the Dirichlet condition but a similar argument is 
applicable to the Neumann boundary condition, corresponding to acoustically hard 
obstacles. 

Fix e > 0, an arbitrary small number. 



4 



SEMION GUTMAN AND ALEXANDER G. RAMM 



Lemma 2.1. There exist L — L{e) and eg. — cg{t) sueh that 

L{e) 

(2.2) IK + 5IcKe)V^||i2(s) <e. 

i=0 

If (2.2) and the boundary eondition (l-l) hold, then 

L(e) 

(2.3) ||we - v\\l^(s) < e, := ^ ceXe)il)e, 

where v is the scattered field defined below formula (1-2). 
Lemma 2.2. If (2.3) holds then 

(2.4) \\\v^ - v\\\ = 0{e) , e^O, 

where \\ \ ■ \ \\ := \\ ■ \\h^ (^o') + II • llL2(_D':(i+|a;|)-T), ^ > I, m > is an arbitrary 
integer, H™ is the Sobolev space, v is the scattered field, is defined in (2.3), and 
both Vc and v in (2.4) are functions defined in D' . 
In particular, (2.4) implies 

(2.5) \\v,-v\\L2iSu)=0{e), e^O, 
where Sr is the sphere centered at the origin with radius R. 

Remark. The proof of Lemma 2.2 in [7] shows that estimate (2.4) for the scattered 
field V in the exterior D' of the obstacle D follows from the boundary estimate 

ll^^e - ^^||z,2(s) < e 

for any outgoing solution of (V^ + k^) = in £)', regardless of how such 

solution f e was constructed. 

Lemma 2.3. One has: 

(2.6) c^(e) -> A,(a), V£, e ^ 0. 

The modified Rayleigh conjecture (MRC) is formulated as a theorem, which 
follows from the above three lemmas: 

Theorem 2.4. For an arbitrary small e > there exist L{e) and c^(e), < ^ < 
L{e), such that (2.2), (2.4) and (2.6) hold. 

See [7] for a proof of the above statements. 

The difference between RC and MRC is: (2.3) does not hold if one replaces Vf, 
by X^fco ^^('^)V'^' ^'^^ L — *■ 00 (instead of letting e — »■ 0). Indeed, the series 
Til 1=0 t diverges at some points of the boimdary for many obstacles. Note 

also that the coefficients in (2.3) depend on e, so (2.3) is not a partial sum of a 
series. 

For the Neumann boundary condition one minimizes 



ON 

with respect to q. Analogs of Lemmas 2.1-2.3 are valid and their proofs are essen- 
tially the same. 

See [10] for an extension of these results to scattering by periodic structures. 
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3. Numerical Experiments. 

In this section we describe numerical results obtained by the Random Multi- 
point MRC method for 2D and 3D obstacles. We also compare the 2D results to 
the ones obtained by our earlier method introduced in [8]. The method that we 
used previously can be described as a Multi-point MRC. It was not an iterative 
method, and its applicability was limited by the problem size, i.e. the number J of 
the interior points xj. It required extensive computational resources (run time and 
memory) needed for the minimization part of the algorithm even for a moderate 
number J. 

The Multi-point MRC method is just the first iteration of the Random method. 
It can be run efficiently with a relatively modest J, where J is kept constant across 
iterations. Also, in the Midti-point MRC method the interior points Xj, j = 
1,2,..., J, were chosen deterministically by an ad hoc method according to the ge- 
ometry of the obstacle D. 

The Random Multi-point MRC eliminates the need for this special procedure. 
See the next section for an additional discussion of this issue. 

As we mentioned previously, [8] contains a favorable comparison of the Multi- 
point MRC method with the Boundary Integral Equation Method, despite the fact 
that the numerical implementation of the MRC method in [8] is considerably less 
efficient than the one presented in this paper. 

A numerical implementation of the Random Multi-point MRC method follows 
the same outline as for the Multi-point MRC, which was described in [8]. Of course, 
in a 2D case, instead of (1.4) one has 

Mx,xj) = Hl^\k\x-Xj\)e'^'^^, 

where {x — Xj)/\x — Xj\ = e*^^. 

For a numerical implementation choose M nodes {tm} on the surface S of the 
obstacle D. After the interior points Xj, j = 1, 2, J are chosen, form N vectors 

n = 1, 2, . . . , iV of length M. Note that N = {2L + 1) J for a 2D case, and N = 
{L -\- 1)^ J for a 3D case. It is convenient to normalize the norm in by 

1 ^' 

Then lluoll = 1. 

Now let b = {ff(im)}m=i, in the Random Multi-point MRC (see section 1), and 
minimize 

(3.1) $(c) = ||b + Ac||, 

for c e C^, where A is the matrix containing vectors a^"^ n = 1, 2, . . . , A'' as its 
columns. 

We used the Singular Value Decomposition (SVD) method (see e.g. [4]) to 
minimize (3.1). Small singular values s„ < Wmin of the matrix A are used to 
identify and delete linearly dependent or almost linearly dependent combinations 
of vectors a^"-*. This spectral cut-off makes the minimization process stable, see the 
details in [8]. 
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Table 1. Normalized residuals attained in the numerical experi- 
ments for 2D obstacles, ||uo|| = 1- 
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Let r"**" be the residual, i.e. the minimal value of $(c) attained after N^ax 
iterations of the Random Multi-point MRC method (or when it is stopped). For a 
comparison, let rJJ™ be the residual obtained in [8] by an earlier method. 

We conducted 2D numerical experiments for four obstacles: two ellipses of dif- 
ferent eccentricity, a kite, and a triangle. The M=720 nodes tm were uniformly 
distributed on the interval [0, 27r], used to parametrize the boiuidary S. Each case 
was tested for wave numbers k = 1.0 and k = 5.0. Each obstacle was subjected to 
incident waves corresponding to a = (1.0, 0.0) and a = (0.0, 1.0). 

The results for the Random Multi-point MRC with J = 1 arc shown in Table 
1, in the last column r™*". In every experiment the target residual e = 0.0001 was 
obtained in under 6000 iterations, in about 2 minutes run time on a 2.8 MHz PC. 

In [8], we conducted numerical experiments for the same four 2D obstacles by 
a Multi-point MRC, as described in the beginning of this section. The interior 
points Xj were chosen differently in each experiment. Their choice is indicated in 
the description of each 2D experiment. The column J shows the number of these 
interior points. Values L = 5 and M = 720 were used in all the experiments. These 
results are shown in Table 1, column t^^. 

Thus, the Random Multi-point MRC method achieved a significant improvement 
over the earlier Multi-point MRC. 

Experiment 2D-I. The boundary S is an ellipse described by 

(3.2) r(t) = (2.0cosi, sini), < t < 27r . 

The Multi-point MRC used J = 4 interior points Xj = 0.7r{^^^^^), j = 1, . . . , 4. 
The run time was 2 seconds. 

Experiment 2D-II. The kite-shaped boundary S (see [2], Section 3.5) is de- 
scribed by 

(3.3) r(t) = (-0.65 -I- cost -I- 0.65 cos 2t, 1.5 sin t), 0<t<27r. 
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Table 2. Normalized residuals attained in the numerical experi- 
ments for 3D obstacles, ||uo|| = 1- 
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III 


1.0 
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3600 
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0.001 


3000 
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5.0 
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5.0 
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5000 
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The Multi-point MRC used J = 16 interior points Xj = 0.9r(^!i^), j = 1, . . . , 16. 
The run time was 33 seconds. 

Experiment 2D-III. The boundary S is the triangle with vertices at (—1.0, 0.0) 
and (1.0, ±1.0). The Multi-point MRC used the interior points Xj = 0.9r(^!^^^), 
j = 1, . . . , 16. The run time was about 30 seconds. 

Experiment 2D-IV. The boundary S is an ellipse described by 

(3.4) r(t) = (O.lcosf, sinf), < t < 27r . 

The Multi-point MRC used J = 32 interior points xj = 0.95r(2ligii), j = 
1, . . . , 32. The run time was about 140 seconds. 

The 3D numerical experiments were conducted for 3 obstacles: a sphere, a cube, 
and an ellipsoid. We used the Random Multi-point MRC with L = 0, Wmin = 
10""'^^, and J = 80. The number M of the points on the boimdary S is indicated in 
the description of the obstacles. The scattered field for each obstacle was computed 
for two incoming directions aj = (^, 0), i= 1,2, where (p was the polar angle. The 
first unit vector ai is denoted by (1) in Table 2, a.i = (0.0,7r/2). The second one is 
denoted by (2), a2 = (7i'/2, 7r/4). A typical number of iterations A^iter and the run 
time on a 2.8 MHz PC are also shown in Table 2. For example, in experiment I 
with k = 5.0 it took about 700 iterations of the Random Multi-point MRC method 
to achieve the target residual r""" = 0.001 in 7 minutes. 

Experiment 3D-I. The boundary S is the sphere of radius 1, with M = 450. 

Experiment 3D-II. The boundary S is the surface of the cube [—1, 1]^ with 
M = 1350. 

Experiment 3D-III. The boundary S is the surface of the ellipsoid x^/16 -|- 

y^ + z^ = 1 with M = 450. 

In the last experiment the run time could be reduced by taking a smaller value 
for J. For example, the choice of J = 8 reduced the running time to about 6-10 
minutes. 

Numerical experiments show that the minimization results depend on the choice 
of such parameters as J, Wmin, and L. 
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4. Discussion of the results. 

Let D be an obstacle, S be its boundary, and uq be the incident field. It is 
proved in [7] that if Ve is an outgoing solution of the Helmholtz equation in the 
exterior domain D' and uq + Vc approximates zero in L^(S')— norm on the boundary 
S, then Ve approximates the exact scattered field v in D', see Lemma 2.2 and the 
Remark after it. The Modified Rayleigh Conjecture approach to obstacle scattering 
problems is based on the following observation: the functions ipe{x,z), z € D 
and their linear combinations are outgoing solutions to the Helmholtz equation 
in the exterior domain. Therefore, one just needs to find a combination of such 
functions that gives the best fit to —uq on the boundary S. Then this combination 
approximates the scattered field everywhere in the exterior D' of the obstacle D 
and the error of this approximation is given in Theorem 2.4. 

Various implementations of the MRC method provide different algorithms for 
the construction of such best fit. The original theory of MRC method given in [7] 
guarantees that one can use all i'eix, z) at z = to obtain the required fit. However, 
it is not always possible to obtain the desired accuracy numerically, because this 
would require a very high accuracy of computations due to the fact that the Hankel 
functions differ from each other by many orders of magnitude if £ is large. Thus, 
one wants to restrict the order L. We found numerically that L = 5 is a reasonable 
value from the numerical point of view. 

To keep L reasonably small, we suggcstcid in [8] to use a batch of interior points 
Xj £ D, j — 1,2 . . . , J, and the associated source functions ipe{x, xj) to find the 
best fit (Multi-point MRC). The results are presented in [8]. They improve with the 
increase in the number J of the interior points. However, for J > 30 the problem 
size becomes too large for our computer system to handle it efficiently, thus limiting 
the applicability of the Multi-point MRC. Still, if one can achieve a satisfactory fit 
with this method (2D or 3D), then it is more efficient than BIEM. Another issue 
in the Multi-point MRC is the placement of the interior points xj. In [8] we used 
an ad hoc procedure, by spreading the points just behind the boundary S of the 
obstacle. 

In the present paper we made a further advance in the MRC algorithm. Our 
Random Multi-point MRC is an iterative method. Thus it allows to keep the 
number J of the interior point relatively small, but still achieve a good boundary 
fit. Undoubtedly, there are many strategies as to how to place the interior points 
Xj. Clearly, they should be chosen distinctly in subsequent iterations. Numerical 
experiments showed that restricting these points to a subset of the obstacle D did 
not produce satisfactory relults. We tried to place the points in the interior of D 
in a random fashion, to assure that the entire interior is utilized. The results of 
this algorithm arc presented in this paper. They show a significant improvement in 
the obtained fit as compared with our earlier method. Also, we were able to obtain 
accurate solutions for 3D problems. Clearly, the random choice of the points is not 
essential to the algorithm. We think that a suitable deterministic choice of points 
could be as successful and are working on finding such a choice. 

It may happen that subsequent iterations bring only a negligible improvement 
to the already obtained minimization values of in the method we proposed in 
this paper, and we do not have a theoretically justified way around this difficulty. 

We are trying to find an optimal placement for the interior points and we hope 
to report the results of this research upon its completion. 
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5. Conclusions. 

For a 2D, or 3D obstacle, Rayleigh conjectured that the acoustic field u in the 
exterior of the obstacle is given by 

oo 

(5.1) u(a;,a) =e^'="-"+^A^(a)V€, i^t := Ye{a')hi{kr), a' = -. 

While this conjecture (RC) is false for many obstacles, it has been modified to 
obtain a representation for the solution of (1.1)-(1.2) and to obtain its error. This 
representation (Theorem 2.4) is called the Modified Rayleigh Conjecture (MRC). 
In fact, it is not a conjecture, but a theorem. 

We propose here an implementation of the MRC method which gives an efficient 
approach to solving obstacle scattering problems in 3D problems with complicated 
geometries. Our implementation of the MRC method worked in the cases considered 
more efficiently than the BIEM method. 

The implementation of the MRC method presented in this paper, the Random 
Multi-point MRC, has been successfully applied to various 2D and 3D obstacle 
scattering problems. This algorithm is a significant improvement over the earlier 
implementation of MRC method, given in [8]. The improvement is achieved by 
allowing the required minimizations to be done itcratively. However, even the 
earlier implementation (see [8]) compared favorably to the BIEM. 

The Random Multi-point MRC has an additional attractive feature: it can easily 
treat obstacles with complicated geometry (e.g. edges and corners). Unlike the 
BIEM, the Random Multi-point MRC can be relatively easily used for solving 
obstacle scattering problems with complicated geometries and rough boundaries. 

Further research on MRC algorithms is under way. The authors hope that a 
numerical implementation of the MRC method will yield a more efficient and eco- 
nomical algorithm for solving obstacle scattering problems than the currently used 
methods. 
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